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Abstract 

We present a quantitative study of the I'N cross section in the neutrino energy range 
10^ < Ey < 10^^ GeV within two transversal QGD approaches: NLO DGLAP evolution 
using different sets of PDFs and BK small-x evolution with running coupling and kinematical 
corrections. We show that the non-linear effects embodied in the BK equation yield a slower 
raise in the cross section for Ey > 10® GeV than the usual DGLAP based calculation. Finally, 
we translate this theoretical uncertainty into upper bounds for the ultra-high-energy neutrino 
flux for different experiments. 


1 Introduction 

The observation of neutrino events with energies in the order of PeV by the IceCube Observatory [1] 
has opened a new era in Neutrino Physics. Currently ongoing and proposed experimental programs 
aim at the determination of the neutrino flux for energies of the incident neutrinos several orders of 
magnitude higher than those measurable at IceCube. Thus, the Auger collaboration has recently 
reported on the limits of neutrino fluxes at ultra high-energy (UHE) in [2]. Also, the LUNASKA 
project [3] performs related measurements using the Square Kilometre Array [4]. Among the 
several astrophysical scenarios proposed as source of UHE neutrinos, the interaction of UHE cosmic 
rays with the cosmic microwave background would result in a non-negligible flux of neutrinos 
impinging on Earth up to energies of Ey ~ 10^2-14 QeV. 

An essential ingredient for the determination of neutrino fluxes from experimental data is the 
precise knowledge of the neutrino-nucleon {vN) cross section, as it relates the neutrino flux with 
the number of observed events at a given neutrino energy. However, the study of the vN cross 
section {(Tyjsi) at UHE implies probing QCD in a kinematic regime unexplored so far in ground 
based accelerators, like HERA or the LHC. This region is characterized by very small values of 
Bjorken-a; [x < 10“^ for Ey > 10^^ GeV) and virtualities of the order of the electroweak boson 
mass squared, Ml ~ 10^ GeV^ (for a more precise discussion of the kinematics, see 
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Section 2). Thus, a reliable calculation of (t^n in the previously uncharted kinematic territory 
relevant for UHE neutrino interactions requires theoretically controlled extrapolations from the 
well tested region of phase space studied in collider experiments. 

Calculating the cross section of any high-energy hadronic collision, like the uN studied in 
this work, requires a priori knowledge of the partonic structure of the proton at all relevant 
observation scales. In practice such information is provided by phenomenological parton hts to 
previously existing data based on the use of perturbative QCD renormalization group equations 
in the framework of well dehned factorization theorems. 

The different QCD approaches for the description of the scale dependence of parton distri¬ 
bution functions - analogously, for gauge invariant operators encoding the parton flux into the 
collision - share a similar strategy of resumming radiative terms enhanced by large logarithms 
to all orders. The most widely used framework are the DGLAP equations [5-7], that is, the 
renormalization group equations that describe the scale dependence of parton distribution func¬ 
tions through a resummation of large logarithms ~ aslnQ^/Qg, with Qo some initial scale. The 
DGLAP equations have been successfully and intensively tested against experimental data and, 
together with asymptotic freedom and factorization theorems, provide a fundamental tool for es¬ 
tablishing controlled theoretical predictions. Successful as they are, the DGLAP equations are also 
expected to break down in some kinematic regimes. In particular, at small values of Bjorken-x, 
large logarithms ~ asln(xo/x) emerge and also need to be resummed to all orders. 

In turn, analogous resummation schemes aimed at describing the small-x evolution of hadron 
structure have also been developed. In this direction in the kinematic (x, (5^)-plane, orthogonal 
to DGLAP evolution, the relevant logarithms are ~ a 5 ln(xo/x), resummed to all orders in the 
BFKL approach [8,9]. Additionally, the enhancement of gluon emission at small-x naturally leads 
to the presence of large gluon densities in the proton and to the need of non-linear recombination 
terms in order to stabilize the diffusion towards the infrared characteristic of BFKL evolution. 
Most importantly, the presence of non-linear terms is ultimately related to the preservation of 
unitarity of the theory. Both the resummation of small-x logarithms and the inclusion of non-linear 
density dependent corrections are consistently accounted for by the Balitsky-Kovchegov (BK) 
equation [10,11], which corresponds to the large-A^^c limit of the B-JIMWLK [11-14] hierarchy of 
coupled non-linear evolution equations. The presence of non-linear terms in the small-x evolution 
equations limits the growth rate of gluon number densities for modes of transverse momentum 
smaller than the saturation scale Qs (see Section 2.2 for a precise dehnition). This novel, semi- 
hard dynamical scale marks the onset of non-linear corrections in QGD evolution and leads to 
distinctive dynamical effects such as the generation of geometric scaling [15]. The ability of the 
BK equation at running coupling accuracy to describe the Bjorken-x dependence of the structure 
functions in Deeply Inelastic electron-proton Scattering (DIS) measured at HERA has been well 
established in a series of recent works [16-18]. Further, it was shown in [19] that BK-based hts to 
low-x HERA data are more stable than analogous DGLAP hts under the change in the boundary 
conditions. The unintegrated gluon distributions resulting from these BK global hts to HERA 
data have found many successful phenomenological applications in the analysis of data from the 
proton-proton, proton-nucleus and nucleus-nucleus experimental programs at RHIG and the LHG 
(see e.g. [20] for a review), providing strong evidence of the presence of non-linear saturation 
ehects in available experimental data. 

The natural question arises of which of these two orthogonal QGD approaches, DGLAP or 
BK, is better suited to extrapolate our well stablished knowledge of parton structure to the 
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region of phase space proben in UHE neutrino of very small-a; and high values. Clearly, the 
reliability of either approach in this intermediate kinematic region cannot be determined on a 
priori theoretical arguments. This is so because, at a parametric level, one expects that both 
large logarithmic corrections ~ InQ^ and Inl/x resumed in either approach to be relevant in 
that kinematic regime. It is also clear that claims in favor of one particular approach should 
not be done solely on the basis of agreement with previous experimental data: it is well known 
that one can obtain an excellent £t to the HERA low-x data with a very reduced number of free 
parameters. Such is a necessary but not sufficient condition. Further, beyond describing existing 
data, the usefulness of a given approach rests on its predictive power towards kinematic regions 
experimentally unexplored so far. This latter condition lessens the reliability of phenomenological 
models not equipped with a well defined QCD dynamical input. In our view the use of one or 
another approach should be considered as a systematic uncertainty associated to the theoretical 
estimate to the uN cross section. It is the goal of this work to provide a precise quantitative 
reference for such systematic uncertainty. As we shall explain in detail in the following sections, 
we find that the differences arising in the calculation of the neutrino-nucleon cross section at UHE 
due to the choice of either the DGLAP or BK approaches are sizeable for energies of the incident 
neutrino > 10® GeV. This differences become as large as a factor 4.5 for the highest neutrino 
energies of Ei, = 10^^ GeV studied here. In line with our expectations, the values of the uN total 
cross sections obtained within the BK evolution approach are systematically smaller than those 
obtained within the DGLAP approach due to the presence of non-linear recombination effects 
accounted for the BK approach. 

This paper is organized as follows. In the next section a brief review of the underlying formalism 
of neutrino deep inelastic scattering is presented. A systematic study of the DGLAP approach at 
NLO has been performed in Section 2.1 including the corresponding error bands both for charged 
and neutral current interactions. In 2.2 we introduce state-of-the art saturation effects in the 
computation of a^N through numerical solutions of the BK equation including running coupling 
and kinematical corrections. Finally, as a phenomenological application we recalculate the limits 
of the UHE neutrino flux obtained by various experiments with our new parametrization for the 
cross section. 


2 Neutrino-nucleon cross section 


The inelastic interaction of neutrinos with nucleons is described in terms of charged current (GG) 
and neutral current (NG) interactions, which proceed through and exchanges, respectively. 
The expression of the differential cross section in the fixed-target frame is [21]: 


Glsf Mf Y 
dxdy vr \Mf + ) 


X 


2t7CC,NC/ r^2\ , M M7CC,NC/ , A 

xy ’ (x, Q ) + (1 - y)E^ [x,Q)+xy (1 
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-^3 



( 2 . 1 ) 


where is the Fermi coupling constant, E,^ is the neutrino energy, s = 2MjqEi, with Mjv the 
nucleon mass and Mi denotes the mass of the charged or neutral gauge boson exchanged. The 
kinematics of this process is described in terms of the virtuality of the gauge boson the 
fraction of the nucleon’s momentum carried by the struck quark Bjorken-a: and the inelasticity 
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y = Q’^/ (a;s). The extension of onr resnlts to the case of nentrino-nncleus cross section is straight¬ 
forward. Nnclear effects like nnclear shadowing or other destrnctive interference effects can be 
neglected at high virtnalities (see e.g. the nnclear PDF sets in [22]). This allows to consider nnclei 
as an ensemble of A independent nncleons, leading to a trivial scaling of the total cross section 
with the atomic nnmber A. Farther, we will assnme an isoscalar target N = (p-h n)/2 which is a 
good approximation to the typical nentrino detector material. 

The integral of the differential cross section Eq. 2.1 is dominated by valnes of ~ = 

10^ GeV^. Hence the typical x valne probed is Xmin ~ Mf /2 MmE^. For Ei, > 10® GeV this 
translates into Xmin < 10“®, that is, for nltrahigh energy nentrinos we are going to work in the 
region of small-x and valnes of Ml w 

Once the nnderlying formalism has been presented, in the following sections we are going to 
compnte the strncture fnnctions Ei{x,Q‘^) in both the DGLAP and BK frameworks. 


2.1 Improved parton model and DGLAP evolution 

In the QGD improved parton model the strnctnre fnnctions Ei{x,Q‘^) are linear combinations 
of parton distribntion fnnctions. Explicit expressions for Ei{x,Q‘^) both at leading order (LO) 
and next-to-leading order (NLO) are given in Appendix A. The contribntion of top sea qnarks 
is completely negligible in the phase space region of onr stndy [23]. So it is the contribntion 
of Eg ’ in the energy range that we are dealing with, therefore we are not going to consider 
it in onr fnrther nnmerical calculations. We have chosen two sets of PDFs, NNPDF3.0 [24] and 
MSTW08 [25]. Both sets of parton distribution function correspond to DGLAP evolution at NLO 
accuracy. Even though the two PDF sets rely on the same theoretical framework, differences in 
the implementation like, for instance, the choice of initial conditions for the evolution, and other 
details result in slightly different values of the parton distribution functions. 

We have found, in agreement with [23,26-30], that the purely NLO, Os-dependent terms in the 
calculation of the structure functions are negligible (< 5%) compared to the LO terms, provided 
that the corresponding linear combinations of PDFs appearing in the leading order expression of 
the structure functions, Eqs. (A.1-A.6), are evaluated at NLO accuracy. We have computed a^N 
in two different ways: using the leading order expressions for the Ei{x, of Eqs. (A.1-A.6) with 
the PDFs evaluated at NLO and with the NLO expressions for Ei{x, Q'^) of Eqs. (A.7-A.12). The 
former is labeled as (8) PDF^'"*^ in Fig. 1 and is the one commonly used in the literature, where 

tts-dependent terms are disregarded. The latter is tagged as (8) PDF^'"® and corresponds to 
the strict calculation of the NLO cross section. The result of these two calculations are shown in 
Fig. 1, left-top plot, whereas their ratio 


^NLO (g, pj3pNLO 

aLO o pdF^^° 


( 2 . 2 ) 


is shown in the bottom plot. The value of R stays close to unity and almost constant for all the 
neutrino energy range studied here. This result is better understood by exploring the relative 
contribution to the differential neutrino-nucleon cross section of each individual parton flavor 
both at LO and NLO accuracies, shown in the right plot of Fig. 1 for = 10^ GeV^, which 
is the typical Q^-value that dominates in this process. Due to the very strong rise of the gluon 
distribution function at small-a; - observed in HERA data and accounted for in the PDF sets 
used here -, it may be expected that this contribution, which only appears in the calculation 
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of the cross section at strict NLO accuracy, would become the dominant one at sufficiently high 
Ey. However, this is not the case as the smallness of the coupling and the convolution with the 
coefficient functions compensate for the growth in the gluon distribution at very small values of x. 
We conclude that the leading terms of the NLO expressions Eqs. (A.7-A.12) provide an accurate 
prediction of the UHE vN cross section. 



Figure 1: Left: Comparison of ® pp)pNLO (jNlo ^ pj^pNLO -yvith their ratio. 

CX f 

Right; PDFs and —C 2 0 PDF of quarks and gluons as a function of x. 

2 % 


In our previous discussion about the kinematical region explored in the neutrino-nucleon in¬ 
teraction we have indicated that high Ey is equivalent to small values of x. The most recent 
update of the experimental data available of DIS [24] sets the x-limit in ~ 10“® — 10“®. The PDFs 
are essentially unconstrained below this barrier and these uncertainties are fully propagated to 
the neutrino-nucleon cross section. Smaller values of Bjorken-a: are accessible in the very forward 
region of the LHC by the LHCf or TOTEM experimental programs, but those data are associ¬ 
ated to very small values of the virtuality < 1 GeV ^ and, hence, fall outside the domain of 
applicability of the DGLAP approach, valid only for perturbatively large Q^-values. In Fig. 2 
we can observe how the x-value providing the largest contribution to the total neutrino-nucleon 
cross section shifts towards smaller values when increasing the neutrino energy. The uncertainty 
in the PDFs in that region is directly translated into the error bands for the cross section, given 
at 68% C.L. As an example, for the highest neutrino energy we hnd that the dominant x-value is 
~ 10“® and the PDFs have a ~ 20% error. 

This uncertainty, however, is not universal for all the PDFs collaborations even though they 
use an analogous procedure: an initial parametrization for the PDF that is htted to data followed 
by DGLAP evolution. Not only the error bands but also the central value changes from one set 
to another, as is shown in Fig. 3. We show the neutrino-nucleon cross section calculated with the 
MSTW08 PDF set in the left plot of Fig. 3. It can be seen how the error bands increase while rising 
the neutrino energy due to the increasing contribution of small-x values reaching a 30% uncertainty 
for the highest energy. The right plot shows the ratio between the cross section obtained with the 
MSTW08 and NNPDF3.0 PDF sets, both for charged and neutral currents. The prediction stemming 
from the NNPDF3.0 set is higher for Ey > 10® GeV becoming incompatible with the MSTW08 set 
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Figure 2: The neutrino-nucleon cross section da/dx with 68 % C.L. error bands for = 
10^’8’i2’i4 GeV as a function of x calculated with the NNPDF3.0 PDF set. 


within the error band at the highest neutrino’s energy. This illustrates the decreasing predictive 
power of the approaches based on DGLAP when extended in the a:-direction to values smaller 
than those included in their htting data sets. Moreover, the MSTW08 PDF set covers values down 
to only X = 10“®, which sets an upper value of Ey = 10^^ GeV for the energies that can be studied 
within this parametrization. In turn, the NNPDF routines provide access to x-values down to 
X = 10“^^ and beyond, allowing to extend the energy range of our study up to Ey = 10^"^ GeV. 
Furthermore, the error bands of NNPDF3.0 are systematically smaller than the MSTW08 ones. Thus, 
in what follows all the DGLAP based results will be calculated using the NNPDF3.0 PDF set. 


2.2 Dipole model and BK evolution 

In the dipole model of DIS at low-a:, the neutrino-nucleon cross section affords the following physics 
interpretation: long before reaching the target, the exchanged electroweak boson fluctuates into 
a colourless quark-antiquark dipole. Subsequently, the qq dipole scatters off the hadronic target 
via multiple gluon exchanges. The DIS structure functions for neutrino-nucleon scattering in the 
dipole model are given by 


77CC,NC 

-^T,L 




91 

Air 


d^bd^r 



i! 


T,L 




r) 


2 

Mdip{x, r, b) 


(2.3) 


where z is the fraction of longitudinal momentum of the gauge boson carried by the quark, r is 
the transverse separation between the quark and the antiquark and b the impact parameter of 
the dipole-target collision (henceforth boldface notation indicates two-dimensional vectors). The 
expressions of the wave functions '0 t,l for the splitting of the gauge boson, with transverse (T) 
or longitudinal (L) polarizations, into a qq dipole at lowest order are given in Appendix B. All 
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Figure 3: Left: the neutrino-nucleon cross section at NLO with MSTW08 PDFs. Right: comparison 
between the results for obtained with MSTW08 and NNPDF3.0. 


the information about the strong interactions in Eq. (2.3) is encoded in A/dip(x, r, b), the dipole- 
hadron scattering matrix. Under the translational invariant approximation, i.e. assuming that 
the nucleon is homogenous in the transverse plane, it only depends on the absolute value of the 
dipole transverse size, r = |r| such that the dipole-hadron cross section is given by 

aqq{x,r) = 2 j d^bArdip(a:,r,b) = 2aoA/'dip(a;,r) , (2.4) 

where ao is a parameter that Exes the normalisation and is htted to experimental data. Although of 
ultimate non-perturbative origin, the evolution of the dipole scattering amplitude towards smaller 
values of x can be studied perturbatively via the BK equation. It reads 

gA/'(r, x) ^ f ri, rs) [Af{ri, x) + A/'(r 2 , x) - A/'(r, x) - N'{ri,x)N'{r 2 , x)] (2.5) 

am(xo/x) j 

where xq is the initial value of Bjorken-x and K, is the evolution kernel, which plays an analogous 
role to the DGLAP splitting functions. The evolution kernel is now known to NLO accuracy in 
a^lnl/x [31] and also at running coupling accuracy [32-34]. However, NLO BK evolution turns 
unstable for a large class of initial conditions [35]. Therefore, here we shall solve the BK equa¬ 
tion either at running coupling accuracy and also adding double logarithmic corrections recently 
calculated in [36, 37] to the evolution kernel. These two evolution schemes proceed through an 
all-order resummation of just a subset of terms appearing at NLO accuracy, thus rearranging the 
perturbative series. 

Similar to DGLAP evolution, solving the BK equation is an initial value problem, i.e. it is 
well dehned only after initial conditions at the initial evolution scale Xq have been provided for 
all values of the dipole size r. This introduces free parameters to be htted to data. Here we 
shall use the parametrizations of the dipole scattering data obtained in References [16,17,38] after 
htting all available data on the reduced cross sections measured in ep collisions at HERA with 
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Figure 4: The neutrino-nucleon cross section calculated within the dipole model using different 
parametrizations for the dipole cross section, aqq: RUN, DLA and PHENO. 



Figure 5: Neutrino-nucleon cross section calculated with NLO DGLAP with NNPDF3.0 and with 
running coupling BK. 




























































X < 10“^. In these works the evolution kernel was evaluated at running coupling accuracy [16,17], 
yielding a very good description of the data y^/d.o.f ~ 1. We shall refer to this parametrization 
as RUN in what follows. It should be noted that a study of the neutrino-nucleon cross section using 
running coupling BK evolution has been presented in [39]. More recently a new parametrization 
of the dipole amplitude tested against the same data set but also including large double transverse 
logarithmic corrections in the evolution kernel [37] (also referred to as kinematic corrections in [36]) 
has been presented in [38]. Again, the reported £t results were very good with y^/d.o.f ~ 1 in all 
cases. We shall refer to this parametrization as DLA in what follows. The main effect of including 
double logarithmic corrections in the evolution kernel is a further slowdown of the evolution speed 
towards small-x. As a consequence the DLA parametrization yields a systematically smaller value 
of the neutrino-nucleon cross section than the RUN parametrization. 

The maximum values of virtuality included in both the RUN and DLA hts is = 500 GeV^. 
Therefore the initial parametrization are not constrained in the region = 10^ GeV^ relevant for 
our calculations. This is an analogous problem to the extension of DGLAP PDF hts to the region 
of very small-x previously discussed. An important difference, though, is that at high rapidities 
Y = ln(a;o/a;), and on account of its non-linear character, the solutions of the BK equation become 
independent of the initial conditions. In other words, the existence of a fixed point in BK evolution 
reduces the dependence on the £t parameters, resulting in asymptotically universal predictions 
that are only controlled by the dynamical information contained in the evolution kernel. A detailed 
study of this feature of BK evolution, referred to as scaling in the literature can be found in [40,41]. 
We have tested the sensitivity to the initial conditions by reducing the maximum included in 
the HERA data htting set from 500 to 50 GeV^. This implies a change in the initial conditions for 
the evolution. We observe that, for a given evolution scheme (running coupling only or running 
coupling plus DLA corrections), this results in negligible difference in the calculation of the uN 
cross sections. 

The evaluation of the total neutrino nucleon cross section involves integration of the differential 
cross section (Eq. (2.1)) over all values of Bjorken-x, 0 < x < 1, whereas the dipole amplitude 
parametrization used here are only well defined for x < 10“^. Other works in the literature 
circumvent this problem by extending analytically the dipole amplitudes into the region x > 10“^. 
Here, in order to account for the contributions to the cross sections from the region 10“^ < x < 1 we 
will assume NLO DGLAP evolution for this range of x. Furthermore, we shall £x the normalisation 
(To imposing that the BK-cross section coincides with the DGLAP one at the lowest neutrino energy 
we are considering, = 10^ GeV. This choice is well motivated since the (x, Q^)-region proben 
at that energy is well described by the DGLAP formalism, that is, large and moderate x. This 
procedure implies increasing the parameter (Tq of the RUN and DLA parametrizations by a ~ 20%, 
from ~ 25 mb to ~ 30 mb approximately. 

For comparison, we shall also calculate the neutrino-nucleon cross section using the GBW [42] , 
HM [43] and Soyez [44] phenomenological models for the quark-antiquark dipole cross section, 
recently used in [27] for the calculation of the vN cross section at UHE. Other phenomenological 
studies of the neutrino-nucleon cross section in different QGD approaches, including the dipole 
model, can be found in ref. [45]. We shall refer to these models collectively as PHENO in what 
follows. Gontrary to the RUN and DLA parametrizations, these models were htted to the old HERA 
data on structure functions, which have been superseded by the much more precise combined 
analyses performed by the HI and ZEUS collaborations [46]. It is also known that more exclusive 
features of the PHENO models are incompatible with more exclusive data like, e.g. single inclusive 
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hadron production cross sections in pp collisions (see [20] for an extended discussion). Finally, 
the energy dependence of PHENO models is mostly driven by a fit parameter, A, that provides the 
x-dependence of the saturation scale Q‘^{x) ~ x~^. Thus, their extrapolation beyond the region of 
phase space covered by the htted HERA data is not supported by a well defined QCD dynamical 
input. 

We present our results for the uN total cross section as a function of the neutrino energy 
calculated within the dipole model in Fig. 4. We observe that the DLA parametrizations yield a 
systematically smaller cross section than then RUN ones. The differences between the two are of a 
factor < 2 for GeV. This is an expected result, since the main role of DLA corrections 

is to further reduce the phase-space for small-x gluon emission, thus resulting in a slower growth 
of the dipole amplitude with decreasing values of x. Once the evolution kernel for BK is fixed, 
the value of is almost insensitive to the initial conditions for the BK equations. This is seen 
by comparing the predictions for either DLA or RUN parametrizations htted to different subsets 
of HERA data, up to = 50 and 500 GeV^ respectively. Finally, the results obtained from 
the PHENO parametrizations yield an even smaller value of the cross section at high energies. We 
conclude that the main systematic uncertainty in the calculation of within the dipole model 
stems from the precise dynamical input embodied in the evolution kernel, and not from the initial 
conditions for the evolution. 

We summarise the results of this section in Fig. 5. There we compare the values for 
obtained under the DGLAP and BK QGD evolution schemes. The band for the DGLAP based 
calculation originates from the uncertainties of the NNPDF3.0 PDF set used here. In turn, the 
uncertainty band associated to the BK-based results is related to the choice of running coupling 
or running coupling plus DLA corrections evolution schemes for the high-energy extrapolation. We 
see that the DGLAP and BK approaches yield similar results for a^N up to energies ~ 10® GeV. 
For higher neutrino energies, the BK-based approaches lead to a systematically smaller value of 
the cross section than the DGLAP based calculation. That difference between DGLAP and BK 
approaches become as high as a factor ~ 4.5 at the highest neutrino energy E^, = 10^“^ GeV. 


3 Limits on astrophysical neutrino fluxes 


In absence of a unihed theoretical framework for the description of QGD scattering in the full 
kinematic (x, Q^)-plane, the differences in the theoretical prediction of the total neutrino-nucleon 
cross section induce an uncertainty in the determination of the upper bounds of astrophysical 
neutrino fluxes. 

The number of neutrino-induced events, in the form of penetrating air showers, at an array of 
water-Gherenkov surface detectors, like those of the Pierre Auger Observatory, is 

TVe^t = J dE,^Stot{E.) (3.1) 


where the total exposure is calculated given the different flavor weigths Ui in the flux, and the 
z/-nucleon cross sections from 


2=152,3 l ' i,Pi 




(3.2) 
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The exposures are the product of the effective aperture, that is the detectors projected area 
weighted by the detection probability integrated over solid angle, and the range of depths within 
which the shower must originate to trigger the device, integrated over time. They are functions of 
the shower energy Egh = yEy and depend on the detection method and the shower type (hadronic 
or electromagnetic, according to the nature of their hrst interaction). Importantly, they can also 
depend on the neutrino-nucleon cross section. In a recent work [2] , a detailed calculation of Auger’s 
total exposure has been presented from the combination of three channels: Earth-skimming neu¬ 
trinos (ES) and down-going neutrinos from high angles (DGH) and low angles (DGL). The hrst 
one involves Vr GG interactions, and the other two all havors with both GG and NG interactions. 
Assuming the same hux of all havors, 0,^/2 = (p^. = pp., with 


dE„ 


= k E-^ 


(3.3) 


and an upper limit of A"up = 2.39 signal events at 90% G.L., from zero observed and zero back¬ 
ground events, including the uncertainties in the exposures, they set a limit 


^90 < 


N, 


up 




= 6.4 X 10 ® GeV cm ^ sr % ^ 


(3.4) 


This limit applies in the energy interval 10® < E^/GeV < 2.5 x 10^° and was obtained adopting 
the DGLAP i/-nucleon cross sections given in [47]. The relative contributions of the three channels 
to the expected rate are (ES):(DGH);(DGL) ~ 0.84:0.14:0.02. According to our running coupling 
BK calculation with DLA corrections (see Fig. 5) their result may be based on an overestimated 
expectation for the event rate. Neglecting in this case the dependence of the exposures on the 
neutrino-nucleon cross-section in Eqs. 3.1 and 3.2, we obtain that the limit fcgo could be enlarged 
by a factor of ~ 1.5, that is the inverse of the ratio between the BK cross section and that of 
DGLAP at an intermediate value of E^ = 10® GeV. 

Over an energy interval, 1.6 x 10® < Ey/GeV < 3.5 x 10®, IceGube [48] has established an 
upper bound on the same flux. 


fcgo < 8.3 X 10 ® GeV cm ^ sr % ^ 


(3.5) 


that should be rescaled similarly. 

Even higher energies is the domain of neutrino radio-detection experiments, based on the co¬ 
herent radiation emitted by secondary charged particles travelling faster than the phase velocity 
of light in a dense dielectric, radio-transparent medium, such as ice, sand or the lunar regolith 
(Askaryan effect). A few experiments operate with terrestrial ice as target, monitored with an¬ 
tennas embedded in ice (RIGE [49]), suspended from balloons (ANITA-2 [50]) or mounted on a 
satellite (FORTE [51]). They are sensitive to increasing energy intervals: 10® < Ey/GeV < 10^^, 
10® < Ey/GeV < 10^^-® and 10^® < Ey/GeV < 10®'^, respectively. In particular, ANITA-2 sets a 
90% G.L. limit to a pure Ep"^ in the energy interval 10® < Ey/GeV < 10^^'® of 

fcgo < 2 X 10“^ GeV cm“^ sr“® s“® (3.6) 

using the DGLAP prediction for the uN cross sections given in [30] . Although the energy explored 
by ANITA-2 is higher than that in Auger, we obtain a similar enhancement factor of ~ 1.4 at 
10^^ GeV since the event rate scales with according to [50]. 
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There are also several experiments detecting radio pulses from neutrino-initiated particle cas¬ 
cades in the Moon, under the LUNASKA project [3]. They offer the largest potential aperture 
but they are sensitive only to the most energetic neutrinos {Ey > 10^^ GeV), so that the pulse 
is bright enough to be visible from such a long distance. Based on a geometric aperture A{E), a 
90% C.L. model-independent limit to a diffuse isotropic flux, 

^ ^^_ ( 3 . 7 ) 

<^Ey Ey tohsAefi{Ey) 

has been presented in [52], where fobs is the effective observing time. An important source of 
uncertainty is the cross section that translates almost linearly to the aperture and hence to the 
flux limit [3]. The flux limits from these experiments can be up to a factor of 2.5-4.5 larger in the 
energy range from Ey = GeV, following BK-based result for ayjsr presented in this work. 

4 Conclusions 

We have presented a quantitative study of the neutrino-nucleon total cross section in two different 
QGD approaches: NLO DGLAP and running coupling BK evolution. We have used state-of-the 
art parametrizations of the parton distribution functions, in the DGLAP case, and of the dipole- 
nucleon scattering amplitude in the BK case. These parametrizations have been successfully 
tested against the available experimental data from HERA and the LHG. The differences in the 
cross sections arising from the use of one or another approach are due to the different theoretical 
input driving the extrapolation to the previously unexplored kinematic territory of UHE neutrino- 
nucleon collisions. The cross sections obtained in the BK framework are systematically smaller 
than those calculated within the DGLAP framework, the difference between these two approaches 
being as large as a factor 4.5 at the highest neutrino energies studied in this work Ey = 10^^ GeV. 
This is an expected result, since the BK framework includes dynamical recombination effects 
that reduce the growth of the gluon densities, and hence of the total cross section as well, at 
small values of Bjorken-x. These systematic differences in the theoretical calculation of the total 
neutrino-nucleon cross section affect directly the experimental analyses for the determination of 
the neutrino fluxes of astrophysical origin. Neutrino experiments exploring the extremely high 
energy range are the best playground to test our predictions, that start to differ significantly 
from those of DGLAP for Ey > 10® GeV. We have shown that for the energies explored by the 
LUNASKA experiments, Ey ~ 10^"^ GeV, the actual predictions for the upper bounds could be 
enlarged by a factor of up to ~ 4.5. 
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/ 

9 I = Ft - QfS\v?0w 

9r = -Q/sin^^iy 

u 

d 

\ — \ sin^ 6w 
— 7 + 1 sin^ 9w 

-|sin^ 9w 

1 sin^ 9w 


Table 1: Relevant couplings for neutral current interactions with 9w the weak mixing angle. 


A LO and NLO DGLAP expressions for pN scattering 

At leading order (LO) we hnd [21] the following structure functions for charged currents, 

2xF^^{x,Q^) = (A.l) 

F^^{x,Q‘^) = 

d) + 5 — c + 6 


— {u -\- d -\- Vj -\- dj^ s c b 


F^^{x,Q‘^) = ^[u + d 
and for neutral current interactions, 


u 


(A.2) 

(A.3) 


2xF^^{x,Q'^) = F^^{x,Q'^) 


NC/ 


^ ^(2 + (5r)^ + {Sp + iOnf] {udud) 

+ 2 M? + [gif] (s + b) + 2 [(gif + (sS)^] c 
Q^) = 1 [( 92 )'^ - (g'hf + (gif - (gif] {a + d-n-d) 


(A.4) 

(A.5) 

(A.6) 


where the couplings are given in Table 1. To account for NLO corrections one must add a 
convolution of the PDFs with appropriate coefficient functions, including a contribution from 
gluons. The structure functions at NLO for Uf active flavors read [53], 


nCCI 


Fr(a;,Q^) = F| 
F^^{x,Q^) = 


+ 


CTs (Q 
27r 


24 r 


C?®-F“|,.o +A/C?®9 


ILO ■ 2 


+ ^4^ [Cl ® F^\^ + UfCl ® xg] 


F^^{xM^) = F^^\ + 


2 ILO ■ 27r 

OlsiQ"^ 


LO 


27r 




F^M^) 

F^M^) 


= Fi 


NCI 


+ 


pNC, , 

2 iLO + 27r 


Cis{Q^) 
LO ' 27r 

dsm 




n 0 + 53[(9{)'^ + (9p'^]Cf 0 9 




Cl ® ^2^^|lo + 2 5^[(f/{)^ + i9Rf]Cl ® xg 
f 


pNCj , 

2^3 Ilo + 27r 


[^q ^ j;,NC 


[C's '<=>' 3 Ilo 


(A.7) 

(A.8) 

(A.9) 

(A.IO) 

(A.ll) 

(A.12) 
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with 




+ (sifl = {2 M? + (saT + 3 M? + (gi?] } (A,13) 

f 

for Uf = 5, where the PDFs are evaluated at NLO. The coefficient functions are given by [21,54], 




1 + 3 1 

-In 2 ;- 

1 — z 


21-z 


+3 + 2z-{l + z) ln(l - z) - 6{1 - z) + 


where 


Af = -2z, A| = 0, A| = -(1 + ;2) 
and for the gluon initiated processes. 


ci(a = ^ {[(1 - + -1 ta ( 


1-z 


-8z^ + 8z-l + A^ 


where 


and 


Al = -Az{l-z), Af = A| = 0 


asm = 


a.(M|) 


&o — 


(33 - 2N,) 


(A.14) 

(A.15) 

(A.16) 

(A.17) 

(A.18) 

(A.19) 


l + a,(M|)/?„ln(QVM|)’ 12i 

A numerically small disagreement with Ref. [53] was found. The plus prescription is defined as 
usuaT and ® denotes the convolution defined as 

(7/0 PDF = [ —Cf{z)PDF(-) . (A.20) 

Jx ^ V^/ 

B Light cone wave functions for the dipole model 

Explicit expressions for the wave functions squared in the massless quark limit are as follows [55] , 

(B.l) 


mir.z^Q'^) = ^[z^ + {1 - z^K^ier) 

TT^ 

4’P(r,-~.Q^) ' = - .-)"0"A'=(er) 


z, Q^) = ^ [{gif + {gif + {gif + {gif] [z^ + (1 - zmKl{er) 


TT 

2 3 


mf, Qf = ^ [{gif + {gif + {gif + {gif] - ^)^Q^A:o(er) 


27r2 

2 24 


where 


= z{l - z)Q^ 

and Kq^i{x) are the modihed Bessel functions. 

dzf{z)g{z)+= [ dz[f{z) - f{l)]g{z) - f{l) f dzg{z). 

J X J X Jo 


(B.2) 

(B.3) 

(B.4) 

(B.5) 
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